Eigenmode description of Raman scattering in atomic vapors in the presence of 

decoherence 
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A theoretical model describing the Raman scattering process in atomic vapors is constructed. 
The treatment investigates the low-excitation regime suitable for modern experimental applica- 
tions. Despite the incorporated decoherence effects (possibly mode dependent) it allows for a direct 
separation of the time evolution from the spatial degrees of freedom. The impact of noise on the 
temporal properties of the process is examined. The model is applied in two experimentally relevant 
situations of ultra-cold and room-temperature atoms. The spatial eigenmodes of the Stokes photons 
and their coupling to atomic excitations are computed. Similarly, dynamics and the waveform of 
the collective atomic state are derived for quantum memory implementations. 

PACS numbers: 42.50.Nn, 42.50.Ct, 42.50.Ar. 



I. INTRODUCTION 

Collective Raman scattering in multi-atomic ensem- 
bles serves as a potential, basic building block for quan- 
tum information processing and quantum state engi- 
neering [i|. It can be used to generate single spin 
waves and photons |2|, |3| , store them Q and shape their 
wavepackets [5j. Recently, a major improvement in the 
control of the atomic ensembles has been achieved by 
means of multiple longitudinal spin wave modes, in- 
dependent light wavepackets can be stored, reshuffled 
and reshaped ja, Hi or interfered with subsequent light 
pulses [»|. Moreover, transverse degrees of freedom can 
be used to generate spatially multimode squeezed states 
[9||. In this article, we formulate a theoretical model 
of the three-dimensional Stokes Raman scattering in 
atomic ensembles both in the ultra-cold and the room- 
temperature regimes. Despite the apparent complexity 
of the problem, we provide a simple method of obtain- 
ing the independently evolving eigenmodes, which fully 
determine the system's dynamics. 

Typically, the Raman scattering process occurs in an 
ensemble of atoms optically pumped into one of the 
ground states, say |0). Then, the collective excitations 
to another ground state, |1), called "spin waves", can be 
utilized to store the quantum information. These excita- 
tions |0) -^ |1) are induced by the ofF-resonant coupling 
of |0) to higher lying states, denoted by \m) in Fig.[Tl ac- 
companied by the spontaneous emission of Stokes pho- 
tons. The strength of the coupling is controlled by the 
Raman pump laser field, Ep. Furthermore, the read-out 
process can be easily performed by means of the inverse 
Raman scattering, in which the information stored by 
the excitations is passed onto the anti-Stokes photons 
produced [l(i| . Here, we investigate the dynamics of the 
information write-down alone. Nevertheless, our model 
can be easily reformulated to describe the read-out case. 

Our treatment is complementar y t o the description 
introduced in the pioneering works (ll| - [l3j . in which the 
macroscopic excitation regime was investigated. In con- 
trast, we provide an exact solution in the low atomic ex- 
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Figure 1. (Color online) Scheme of a Raman transition. 
An atom in state |0) is illuminated with an external pump 
field Ep, detuned from the |0) — > \m) transition. A single 
collective excitation, quantified by the operator b^ , is created 
by the |0) — >■ |1) atomic transition. It is achieved owing to 
the absorption of a photon from the pumping beam and 
emission of a Stokes photon described by the field operator 
Es. 



citation number limit, especially important for modern 
applications. Notably, we address the typical Gaussian 
transverse pumping beains, as opposed to the flattop 
profiles considered in |ll| - |l3j . Furthermore, rather than 
focusing on the near field, we directly establish the far- 
field properties of scattered light, measurable in current 
experiments. In our approach, we neglect higher-order 
self-interaction effects, which lead to the build-up of the 
excitation gradient along the sample (l3l . [l4| . Thanks 
to this simplification, which is valid in the regime of a 
low number of atomic excitations, we are able to de- 
couple the problem's spatial and temporal degrees of 
freedom. This enables an in-depth analysis of the mul- 
timode capacity of the process. Hence, we provide the 
estimates of its eigenmode parameters, describe their 
evolution and include the possibly mode-dependent de- 
coherence effects. In particular, we can easily account 
for the fact that diffusion most efficiently destroys the 
spin waves characterized by the short-distance spatial 
phase and amplitude variations. Crucially, we describe 
the final atomic state, which, serving as the quantum 
memory, could potentially be mapped onto the light in 



a read-out process (T5l.ll6|. 

This article is organized as follows. In Sec. [IT] we de- 
velop a general three dimensional model to obtain an 
eigenmode description of the system. In the eigenba- 
sis, the Heisenberg equations coupling the atomic exci- 
tations with Stokes photons describe the pairwise cou- 
pling of the atomic and field modes, preserving the inde- 
pendence between distinct mode pairs. Moreover, this 
approach allows for the simple incorporation of the de- 
coherence effects, typically diffusion and spontaneous 
emission. The Raman gain for each spin wave opposes 
the decoherence; thus by varying the pump intensity 
and the size of the sample one can change the number 
of spin-wave modes that are macroscopically occupied. 

Section IIIII discusses the temporal evolution of the 
statistics of each eigenmode. We show that, even for the 
low-occupied spin wave modes, the number of atomic 
excitations can be estimated with a precision that over- 
comes the quantum projection noise by counting the 
scattered photons. Hence, these results may facilitate 
experiments, in which the transverse multimode capac- 
ity of the scattering [9] would be combined with the 
storage 0|- 

Finally, in Sec. IIVI we perform the spatial singular 
value decomposition that is valid in two typical experi- 
mental situations (i.e. when an ultra-cold atomic cloud 
is interacting with a uniform Raman pump or when 
a room-temperature atomic cell is illuminated with a 
Gaussian pumping beam). 



II. THE MODEL 



When N is large, one can replace the summation over 
separate particles in Eq. (|2b[) with an integral over the 
atomic density n{r). Within this approximation, the 
Hamiltonian of the Raman scattering process reads 

H = ih [fdk dK h{k, K) at (k, t) 6^ (K, i)e*'^'=* + H.c, 



(3) 
with ojfc = ck — ujs. where we have introduced the central 
Stokes frequency uJs = ^p — woi- The function /i(k, K) 
determines the probability amplitude of the photon-spin 
wave pair creation with wave vectors k and K, respec- 
tively, and reads 

/i(k,K) = -^^ Jdre-^^^+'^-'^'-y'^y^Spir), (4) 

where go is the coupling constant [l^. Note that the 
pump amplitude £p has been chosen to be time inde- 
pendent, which corresponds to a common physical situ- 
ation of square pulses. As this assumption is made only 
to simplify the calculations, one could in principle gen- 
eralize the model for cases, when the pump amplitude 
£p{r, t) is a separable function of r and t. 

The Hamiltonian ^, together with the commutation 
relations, gives the set of coupled Heisenberg equations, 
describing an evolution of the multimode fields of pho- 
tons and atomic spin waves, 

dtd{k,t)^ fdKe"^''*h{k,K)b\K,t) (5a) 



(5b) 



dtb{K,t)^ /dke*'^'=*/i(k,K)aT(k,t) 



We consider an ensemble of identical atoms, each be- 
ing initially in its ground state |0). Owing to the in- 
teraction of the probe with an intense, non-resonant, 
linearly polarized Raman pump beam with amplitude 
£p, frequency ujp and wave vector kp, that is, 



^j,(r,i)=£p(r,i)e^(''''-"-"''*)-t-i 



(1) 



some atoms are transferred to the other ground state 1 1) . 
Each atomic transition is accompanied by an emission 
of a Stokes photon, as schematically depicted in Fig. [TJ 
To analyze this process in detail, we describe the 
quantum state of Stokes photons and atoms respectively 
using two field operators 



£:(+)(r,i)= /dke'(''"-'='=*)a(k,i) 



b(k,t) 



7^^' 



,i(k-rQ— cjoit) 



(2a) 

|0)a(l|a, (2b) 



where c is the speed of light, luqi is the |0) — ?► |1) transi- 
tion frequency and k = |k|. The operator b annihilates a 
quantum of the collective atomic excitation (i.e. a spin 
wave with momentum hk) . The index a runs over all the 
N atoms in the ensemble. If almost every atom occu pies 
the |0) state, the Holstein-Primakoff approximation jlTJ 
holds, so b satisfies the bosonic commutation relations. 



A. Modal decomposition 

The modal decomposition of the fields is achieved by 
an adequate change of the basis for photons and spin 
waves. To demonstrate this, we start by integrating 
Eqs. ((5a|) and (j5b|) up to the first order in the coupling 
constant go- The integral of exp(icL'fct)/i(k, K) over an 
interval t G [0,r] is proportional to sinc(^^) /i(k, K). 
Let us note, that the sole h{k, K) corresponds to an in- 
finitely short interaction time (r— >-0) and, in principle, 
would give rise to scattered photons of arbitrary fre- 
quency. Moreover, as indicated by Eq. Q, the width Wk 
of h{k, K) as a function of k is inversely proportional to 
the effective spatial spread of the product ^/rUr) £p{r) . 
Hence, as in typical experiments the atomic sample is 
less than 10-cm long, Wk is constrained to be greater 
than 10 m^^. On the other hand, the sinc(^^Y^) part is 
narrowly peaked around ks =uJs/c, so that even for an 
extremely short interaction time of r '^ 1 ns its width 
(ct)^^ ^ Wk- Thus, for any real interaction time, it is 
the sine function that defines a thin spherical shell of al- 
lowed k vectors for Stokes photons with radius k k, kg. 
Therefore, in spherical coordinates, the radius of the 
first argument of h can be fixed at kg. In other words, 
since the interaction lasts longer than the time it takes 



for photons to travel through the sample, the spectral 
width of the scattered light is independent of the shape 
of the atomic cloud. This is a significant step, which 
allows one to perform the singular value decomposition 
|19| as follows: 



smc^- 2 



/i(k = {kKks,k^,ke 



,K) 



* Irn 

with the k dependence remaining only in the spectral 
factor. The integer numbers [Im) parametrize the dis- 
tinct modes and pair the mode functions of photons, 
?/'J'^(fc^, /cg), and spin waves, V'fm(K), which both form 
separate orthonormal bases. The real and positive sin- 
gular values, denoted by C,im, quantify the coupling 
strength between V'J^ and ipf^ for each {Im). The de- 
composition presented in Eq. ^ is always possible for 
any regular coupling function and does not depend on 
the time interval r provided it is sufficiently long. 



B. Uncoupled mode basis 

Using the mode functions from Eq. (jB]), we expand the 
Stokes photon and the atomic spin wave annihilation 
operators respectively as 



a(k, = X! '^^nSK^^s) O'hnik, t) (7a) 

(7b) 



Irn 



6(K,i)=^^t(K)6,„(i), 

Ira 



where we implicitly assume that the length of the pho- 
ton k vector varies only slightly around kg. These ex- 
pressions substituted into Eqs. (|5al) and (|5b[) give a pair 
of coupled equations for each independent mode {Ini), 
that is. 



dtaim{Kt) ^ Clrae^^^'bln^t)^ 



(8a) 



dthm{t) = Qimjdke-'^'^'dl^ik.t). (8b) 

Note that by taking the Fourier transform of Eq. ([5a|l in 
fc, it becomes apparent that the interaction of the spa- 
tially extended Stokes field with the atomic excitations 
is modeled to be point-like. Furthermore, any higher or- 
der effects, such as the possibility of atomic excitations 
being carried over from one part of the sample into an- 
other, are excluded. In return, since the spin wave mode 
functions are time independent, we manage to separate 
the spatial and temporal degrees of freedom. Precise nu- 
merical investigations in one dimension [lj| show that 
our approach is justified, as long as only a few excita- 
tions are present in the sample. 



C. Input-output relations of a mode pair 

The Eqs. (|5a|) and (|8b[) resemble the analytically 
solvable Wigner-Weisskopf model of spontaneous emis- 



sion. Their solution in time can be derived employing 
the "time slicing" method explicitly explained in Ap- 
pendix |X1 By introducing the operators denoting the 
Stokes field entering and leaving the atomic sample at 
time t, aY^{t), and a°^{t), the fc-dependence is fully 
dropped and the final result reads 



&Z\t) 



dr 



KAt) 



S{t -t)+ CLe 
+ Cin.ei'^--'bljO) (9a) 

Jo 

+ e^^-"*&i„(0). (9b) 



The above set of equations, which preserve bosonic com- 
mutation relations, describes the process of the coupled 
amplification of Stokes and the atomic fields. The pho- 
tons and spin waves are created in pairs with a rate 
defined by the singular values Qm of the coupling func- 
tion h, see Eq. (|6]). This statement is no longer true 
in the presence of losses, which impact on the system's 
dynamics we investigate in the following section. 



D. Atomic decoherence 

The decoherence, which affects the atomic ensem- 
ble, is normally described by a combination of damp- 
ing, caused by the random thermal motion and spon- 
taneous emission of atoms, and space-uniform (illumi- 
nation independent) losses, which are, for example, a 
consequence of atomic collisions. Neither the random 
thermal motion nor the uniform losses mix the distinct 
(Im) modes. The rate of spontaneous emission, how- 
ever, which is a function of the pumping beam's inten- 
sity, might vary between various parts of the sample. 
This, in principle, may lead to a coherent deformation 
of an atomic mode (i.e. a coherent cross-talk between 
different spin waves). In our analysis, we assume such 
effect to be negligible. 

The decoherence process, which does not mix 
the modes, is modeled by adding an expression 
-TimUmit) + V^r^fimit) to Eq. ^. The first term 
represents an exponential decay at the rate Tim, whereas 
the second part is the Langevin noise and ensures the 
conservation of the commutation relations of both the 
photonic and atomic operators, since [/;,„. (t), ///„j/ {t')] = 
S(t-t')dwS„i„i'- 

The updated Eqs. ([5a|) and (j8bl) can be integrated, 
giving the input-output relations of the Raman scatter- 



ing in the presence of atomic decoherence, that is, 



N 



ZHt)^ dr S{t-T) + CLe 



2 ^7i™(t-r) 



atir) (10a) 



+ Ome^-* bjjO) + 0„,V^^ / dT e-^-(*-^)//jr) 



&i™W=0™/dre^-(*-^)[ai^(r)]^ 



(10b) 



+ e^""*6,„,(0) + ^21^^/ rf^e^'"^*""^//m(^), 

where 7/™ = Cf^/2 - Ti^. 

We expect a threshold behavior: When Cf™ > SF;™, 
the pair creation rate exceeds the losses, giving exponen- 
tial growth in time of the population of the atom-spin 
wave pairs. In the next section, based on Eqs. ()10a|) and 
(jlObp . we discuss in more detail both the temporal and 
spatial coherence of the system. 

Finally, let us briefly comment on how our result can 
be extended to describe other types of Raman scatter- 
ing. For example, if the pump's frequency was tuned 
to level |1) instead, the interaction would lead to an 
exchange of excitations between the light and atoms, 
as generated by an a^b type Hamiltonian rather than 
the one of Eq. ^. Thus, one can construct the cor- 
rect input-output relations from our result in a straight- 
forward manner by removing the daggers from all the 
operators in Eqs. (|10ap and (|10bp and substituting 
Cim -^ iC,im- Moreover, in the case of the pumping beam 
coupled to both |0) and |1), following a similar proce- 
dure, one would recover the equations described in |.20.] . 



III. TEMPORAL PROPERTIES OF PHOTONS 
AND ATOMS 

Since the Hamiltonian Q is quadratic and the initial 
state is a vacuum, then all the temporal properties of 
the system are fully determined by the following set of 
correlation functions: 



(at(t)a(i'))-^((F + 7)e^(*+*') 
7 V 



C 

7 



Fe"^! 



t-t' 



c 

7 



,7(t+t': 



atlt-t' 



(bit) a{t')) = - ((F + 7)e^(*+*') - Fe^l*-*'!) 



(11a) 
(lib) 
(lie) 



Because the separate modes of Eqs ([TU]) evolve inde- 
pendently, the mode indices (Im) were dropped in the 
above equations without loss of generality. Also, as only 
the output photons are considered, we have introduced 



a{t) 



\t). 



In the decoherence-free case, when F = 0, the three 
correlation functions pT|) factorize. Hence, all the 
temporal properties are determined by a single mode 
ex exp[((^^t)/2]. Moreover, at any time t, the average 
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Figure 2. (Color online) The mean number A^ of scattered 
photons {dotted lines) and spin waves (dashed lines), as a 
function of t, for C = 1 a-nd three different values of decoher- 
ence strength F. When F = 0, the numbers of the atomic 
excitations and Stokes photons are equal {black dash-dotted 
line). The curves flatten with increasing F and for strong 
enough damping, here F > 1/2, the number of spin waves 
saturates {blue dashed line). 



cumulative number of scattered photons is equal to the 
mean number of atomic excitations present, that is. 



{b\t) h{t)) = dT {ciHt) d{T)) = e'^'* - 1. 



(12) 



This expression indicates that the system can be de- 
scribed as a sum of independent squeezed states for each 
{Im) mode, since Stokes photons and spin waves scatter 
in correlated pairs. 

The situation changes dramatically in the presence 
of decoherence, when F ^ 0. The correlation functions 
pT|) no longer factorize, thus the fields a{t) and b{t) 
become multimode. First, we investigate the impact 
of decoherence on the mean number of scattered atoms 
and photons for various values of F with fixed C = 1, 
see Fig. [21 Clearly, with the strength of decoherence 
rising, the rate at which excitations are created drops. 
Furthermore, the mean number of scattered atoms be- 
haves in two distinct manners. When F < C^/2, then 
7 > and {b\t) b{t)) grows exponentially in time. It 
saturates, on the other hand, for F > C^/2, which can 
be easily inferred from Eq. (|llb[) . On the contrary, the 
average number of Stokes photons obtained by integrat- 
ing Eq. piap increases irrespective of F. 

Second, to examine the effect of decoherence on the 
system in more detail, we focus on the normalized tem- 
poral correlation function of photons. 



9i{t,t') 



(at(t)a(t')) 



vmfmjWWTWYJ' 



(13) 



which can be directly accessed in an experiment. In 
the lossless case, only one temporal mode is present, 
hence for any t, t': gi{t,t') = 1, which confirms the 
full time coherence of the system. When decoherence 
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Figure 3. (Color online) The normalized temporal corre- 
lation function of scattered photons, g\{t,t'), for the sub- 
critical and supercritical losses regimes [i.e., T < 1/2 and 

r > 1/2, (c = 1)]. 



effects are included, the two distinct regimes are again 
observed. In the subcritical case of F < C^/2, [e.g., 
Fig.[3l^a)], the losses are not large enough to completely 
suppress the gain. Therefore, gi always increases up to 
1, as one moves upwards along the antidiagonal direc- 
tion on the diagram. On the contrary, in the supercrit- 
ical losses regime of F > ("^ /2 [e.g.. Fig. M.b)], the full 
time coherence only persists on the antidiagonal [i.e., 
for alH > 0: gi (t, t) — 1] and the correlation function al- 
ways exponentially drops in the antidiagonal directions. 
Looking back at Eq. (|llap with 7 < and calculating 
lim(f_|.f/)_j.oo gi(t, i') = e"*"'*"* ', we realize that, in the su- 
percritical regime, gi decays away from the antidiagonal 
even infinitely far from the origin t = t' = 0. 



A. Estimation of the number of atomic excitations 

In a typical scattering experiment, the number of 
Stokes photons is measured as a function of time. Based 
on this information, one can estimate the number of 
atomic excitations created in the process. Here, we 
demonstrate how to maximize the precision of the esti- 
mation procedure by utilizing the correlation functions 

ofEq. in]). 

To construct such an estimation strategy (optimal for 
any instance of an experiment of duration T) we need 
to minimize the mean squared difference between the 
probability distributions describing the real and the es- 
timated numbers of spiir waves created in the process. 
To this end, we construct an operator representing the 
statistical distance between those distributions, that is. 



V^b\T)b{T)- / dta\t)a{t)w{t), (14) 

Jo 

where differeirt estimatioir procedures correspond to var- 
ious time-depeirdent weight functions w{t). Then, the 
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Figure 4. (Color online) The optimal weight function w(t) 
for estimating the number of atomic excitations present at 
time T, when the pumping beam is switched off, ((" = 1). In 
the lossless case of F = 0, «;(i) = 1, as the average numbers of 
spin waves and detected Stokes photons coincide at all times 
and no estimation is required. The optimal weight function 
becomes steeper, as the strength of decoherence rises. 



optimal strategy is represented by w{t), which mini- 
mizes the mean squared error, that is, 

{V'^) = nb{nb + 2) + / dtdt' w{t)w{t'){a\t)a{t')) 




dtw{t)[w{t) - 2nb]{d\t)d{t)) 



2 dt\{b{T)d{t))\^w{t) 
/o 

r \ 2 

dtw{t){d\t)d{t)) 



(15) 



with rib = {^KT) b{T)) being the mean number of spin 
waves present at T. Hence, the optimal w(t) should 
satisfy the differential criterion 5{T)'^)/5w = 0, which 
corresponds to the equatioir 



dt' {\{d\t)d{t'))\^ + {d\t)d{t)){d\t')d{t'))) w{t') 

-wit){a\t)d{t)) = |(6(r)a(t))|2 + nb{d\t)dit)) (16) 



that is fully determined by the correlations functions of 
Eq. (HH). 

The above equation can be easily solved numerically 
by approximating all of the terms piecewise liirearly. 
The result is a rising function, as illustrated iir Fig. U) 
Its shape can be explained by realizing that for a pump 
pulse of duration T we expect some excitations pro- 
duced at the begimring of the interaction i ~ to wash 
out, while those produced at i ~ T to generally prevail. 
Therefore, it is indeed judicious to observe the photon 
counts registered at the beginning less than those pro- 
duced at the end of scattering. 

It is interesting to note that the application of the 
proper weight function will lead to a sub-Poissonian 



(sub-projectioii-noise) precision of the estimation of the 
number of atomic excitations (i.e. (P^) < rih) for any 
setting. For example, although for large losses (^ — 1, 
r = 1, and very long pulse duration T > 10 the mean 
number of atomic excitations and the average total num- 
ber of scattered photons are only vaguely correlated, 
using a weight function ■w{t) — 1.24exp[2.23(t — T)] 
one can estimate the number of atomic excitations with 
noise over 3-dB smaller than the projection noise. 



IV. SPATIAL PROPERTIES OF PHOTONS 
AND ATOMS 

To obtain a complete picture of the Raman scatter- 
ing process, one needs to know the spatial functions 
i^f^ikip, kg) and t/jf^^CK) used in the original (Im) mode 
decomposition of Eq. ([7]). These mode functions, of 
Stokes photons and spin waves respectively, were de- 
fined in Eq. ^ by the singular value decomposition 
of the coupling function /i(k, K). In this section we 
demonstrate how they can be easily found in typical ex- 
perimental conditions. We focus on the two situations: 
The case of cold atoms, when the shape of atomic cloud 
is characterized with a Gaussian distribution, and the 
room-temperature case, when atoms are uniformly dis- 
tributed in a cell. Although the final form of ipfm and 
ipf^ has to be numerically computed, its accuracy is as- 
sured by many steps that are analytically performed in 
the construction procedure. 



A. Coupling function h(k, K) 

The coupling function h is defined in Eq. ^ via a 
Fourier transform of the product •\/n(r) £p{v), where 
n{r) stands for the atomic density and fp(r) is the am- 
plitude of the classical pump. Since the system has a 
rotation symmetry around the direction of the pump- 
ing beam, taken to be z, it is sensible to describe the 
process in cylindrical coordinates r = (p= {p, Lp),z). 

(a) First, for the case of atoms at low temperatures^ 
we assume their density to be distributed according to 
a Gaussian, that is. 



n{v) = n{p, z) = 



nV 



(17) 



where V is the volume of the probe, nV is the to- 
tal number of atoms, and cTa.p\z effectively parametrize 
the spread of the atomic cloud. The amplitude of 
the pumping beam is taken then to be uniform, hence 
Epijc) = i/io/eo with /q being the pump's intensity. 

(b) Second, in the room-temperature case, we assume 
the atomic density to be constant and equal to n, but 
confined to some finite axial region \z\< L/2 with L be- 
ing the length of the atomic cell. We take the amplitude 



of the pumping beam to be Gaussian, hence 






Up 



(18) 



where ^p — 2z/ (kpa'i) and Op is the beam's waist size. 
Both cases lead to the same general coupling function 



/i(k,K) 



.90 



ft(27r)3 



^^g-.(k+K-k,).rg- 



'Z{z), 



(19) 

for which, to recover the regime (a), one has to set: 
^ = and a = (Ja,p, whereas for (b): S, — S,p and a — Up. 
The function Z{z) is then defined as follows 



Z{z) 



nip 



i^y 



1 



for (a) 
for (b) 



(20) 
In real experimental setups, one calibrates the detec- 
tion system to focus on the Stokes photons that are 
most probably produced in the process but still can 
be efficiently measured. Normally, these are the ones 
produced along the atomic sample at small angles to 
the direction of the pumping beam. Hence, using again 
the fact that the length of the Stokes photon wave vec- 
tor is approximately fixed at fc = |k| « fcg, we ap- 
ply the paraxial approximation as kz ~ ks — k'^/{2ks) 
in the reference frame, in which the atomic wave vec- 
tor transforms to K — >■ K -|- kp — fc^e^. Moreover, 
to diminish the number of free parameters, we intro- 
duce the dimensionless variables in the position space: 
z ~ azZ, p = \J azjks p and similarly for the wave vec- 
tors /« = {k, K}: Kp = yJksJGz kp, k^ = (l/az)^^, with 
Gz quantifying the effective size of the atomic ensemble 
in the z direction. 

Finally, we obtain the revised coupling function 



/i k,K = 



50 



;i(27r) 



^f g-4(kp+Kp)-p+(/?.-fc^/2)5] 



^-p^ l(2-Kp\a,a^\s+i2Kz) 



Z{z), 



(21) 



where now, to retrieve the case (a), one has to set: a = 
'^a,p, Oz — <ya,z and K = 0, whereas for (b): a = Up, 
Gz — L and k — kg/kp w 1. The transformed Z{z) 
function now reads 



Z(5) = 



I nip 
epkl 



\[{i)^e-^\ for (a). 



(22) 



i^ 



for (b). 



Most importantly. Eq. (1211) up to a multiplicative factor 
depends only on a single parameter, the Fresnel number, 
of the illuminated portion of the sample 



F = 



kscr'^ _ ks_ j (yl^p/(Ta,z , for (a) 



(23) 



2TTaz 27r \cr;/L, for (b), 

which is fully determined by the geometry of the probe. 



In conclusion, having fixed the value of F and per- 
formed the singular value decomposition of ft,(k, K) de- 
fined in Eq. (j2ip . we can simply establish the required 
modes of /i(k, K), see Eq. ([5]), just by shifting back the 
K vector and rescaling the obtained eigenmodes. 



B. Mode functions of photons and atoms 

We are now in a position to evaluate the spatial mode 
functions ipf^.^ and ipf^^ of the Stokes photons and atomic 
excitations used in Eq. ([7]). First, we express them in 
orthonormal functional bases. In that way, the singular 
value decomposition of the coupling function /i(k, K) 
can be numerically performed on a discrete tensor of the 
basis coefficients. As argued in Appendix [B] this tensor 
can be computed analytically in the position space using 
the Laguerre-Gaussian modes as the basis functions for 
both TJjf^ and i/'frn- ^^ *^^ result of this procedure, we 
obtain the Fourier transform, J^, of Eq. ([5]): 



(a) Atoms at low temperatures. 



T 



^sinc(^) Mk,K)j =^0™^j:::(r) Vr™(r'), 



(24) 
where the mode functions in the paraxial approxima- 
tion, owing to the cylindrical symmetry of the left-hand 
side, can be written as 



ph/'^^ _ 



«(r) 

V'zt(r) 



■^, 



ph 









(25a) 
(25b) 



Their relation to the original ■(/'^^^^(k) and ■0;^^(K) can 
be found in Appendix iBl 

Figure [5] shows the dominant singular values Qim for 
the Fresnel numbers F — {0.1, 10} in both (a) and (b) 
regimes. In the ultra-cold atoms case, (a), we restrict 
ourselves in numerics to the axial range \z\< •iaa,z, so 
that the Gaussian distributed atomic cloud is well inside 
the sampled region. On the other hand, at room tem- 
peratures (b), as we consider Stokes photons scattered 
nearly with the momentum of the pumping beam, we 
set K = ks/kp = 1. Furthermore, in this case we plot the 
first three photonic and the first atomic eigenmode for 
F = 10. In both regimes, we explicitly write the pri- 
mary unnormalized singular values, CoOi with numerical 
precision to the third decimal place, which are related 
to the previously introduced (^oo via rescaling, that is. 



Coo — Ci 




(26) 



with As = (27r)//cs. 

Note that for large Fresnel number cases {F = 10) in 
Fig. El Cim^s decrease slowly with the mode indices (Im) 
increasing. This can be easily understood considering 
the geometric interpretation of F ex cP' ja^. Large val- 
ues of the Fresnel number correspond to moderate axial 
lengths, (Tz, of the effective interacting region. In this 




(b) Atoms at room-temperature, 



1.38S 




Figure 5. (Color online) The singular values (";„ of dominant 
(Im) modes calculated for the Fresnel numbers F = {0.1, 10} 
in the cold atoms and the thermal atomic vapor regimes. 
In the latter case and F = 10, we also plot the photonic 
eigenmodes 4'JrniP^^~^) ^^^ m = 0, I — (solid blue line), 
7n = 0, / = 1 (dashed blue line) and m = 1, I = (dotted 
red line) , and the modulus of the first atomic mode function, 
|^oo(P,^)!. 



case, many modes fit the illuminated portion of the en- 
semble of diameter a. Inverting this argument, for small 
F only the most fundamental modes are strongly cou- 
pled and contribute to the scattering process. Hence, we 
observe a steep drop in the singular values for F = 0.1 
in both regimes. 



V. CONCLUSIONS 

We have developed a theoretical model of three- 
dimensional Raman scattering process in atomic vapors, 
which is applicable when the number of atomic excita- 
tions created is low. Because of its simplicity, it allows 
to immediately separate the system's spatial and tempo- 
ral degrees of freedom. Moreover, with the decoherence 
effects included, we have found the time evolution of 
the numbers of both the Stokes photons and the spin 
waves produced in the process. Furthermore, by in- 
vestigating the temporal coherences of the system, we 
have demonstrated that despite losses the number of 



collective atomic excitations can be estimated from the 
time evolution of the light field with a precision that 
surpasses the projection noise limit. For the cases of 
ultra-cold and room-temperature atoms, we have shown 
how to easily utilize our model to evaluate the spatial 
mode functions of both the Stokes photons and the spin 
waves. Although the final form of the eigenmodes had 
to be computed numerically, we have shown that, when 
either the atomic density within probe is uniform and 
the pump beam is of a Gaussian shape or vice versa, 
the calculation can be analytically pursued till the final 
decomposition by means of working in the basis of the 
Laguerre-Gauss functions. 

Finally, we point out that the above model can be 
easily generalized and converted to describe other types 
of Raman scattering. Most importantly, it can be sim- 
ply modified to represent the inverse read-out process, 
in which the emitted photons carry the information that 
was originally stored by the spin waves already present 
in the atomic ensemble. 



we obtain the field operators: 



-out _ -in 
"'Im.j IrnJ 



-1 +VT^Chni 1 



^Im^ \ Lt 



■'lm,0 



J-2 / .2 ^\J"»-2 



(Afa) 



Hm.'j 



Ulm.O 



+ ^CimJ2[^ + ^) [«^™.^] ■ (Alb) 

i=0 ^ ^ 

Effectively, the whole field of scattered photons is di- 
vided into slices of length ct. Although at each step, 
the solutions are perturbative, we obtain an exponential 
growth of the effective number of the Stokes photons 
produced due to the constant accumulation of atomic 
excitations over time. 

Finally, we construct the field operators, which con- 
tinuously change in time, that is. 
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Appendix A: "Time slicing" method 



To derive a non-perturbative solution of Eqs. (j8al) and 
(j8b[) , we first integrate Eq. (j8b[) perturbatively in an in- 
terval t e [0,r]. The result obtained is substituted into 
Eq. (I8a[) . which is then integrated in the same time in- 
terval. This way, we obtain the output field of scattered 
photons, denoted by a°^-^, together with bim,i that are 
valid at time t. At the beginning of the next interval 
t G [t, 2t] , we assume that another unperturbed Stokes 
field aJJ^ 2 contributes to the process, as if it entered the 

probe. This operator together with bim.i defines the 
new initial conditions for the repetition of the proce- 
dure in the next "time slice": t £ [t, 2t] . Replicating the 
method iteratively, up to the j'-th consecutive interval. 



a;::f°-(t) = hm 



'J -injout 






bim{t) = lim bin 



]-HX 



hence we obtain the set of non-perturbative Eqs. (|9a| 
and (l9bl). 



Appendix B: Singular value decomposition (SVD) 

To discretize the description of the system, we seek 
a functional basis, in which the eigenmodes tpfmi^) '^^'^ 
V'f^j(r), introduced in Eq. (|24|) . could potentially be ex- 
pressed. Investigating, their definition we explicitly per- 
form the Fourier transform of Eq. ^ . Focusing on the 
photonic (Im) eigenmode, we obtain 

Jo Jo 

(Bl) 



X 



J2Km'ike)e^'^"^', 



where we expanded ipf^ in a multipole series and, as be- 
fore, assumed that the spectral factor just fixes fc= |k| « 
kg. Furthermore, we restrict ourselves to Stokes photons 
produced in the positive direction at small angles, for 
which ipfjn-m' ^^ narrowly peaked around kg = 0. There- 
fore, we can apply the paraxial approximation, in which 

sin kg Ki kg, cos /eg sa 1 — # , and extend the polar angle 



integral's range to kg G [0, oo[. Hence, with a^tp—k, 



v 



/,ph 



E' 







eKge 



27r 



/ 7 ik^kfiCQsap iam' 

X dae ' e 



«(r) 



where according to the Hankel transform 



/,Ph 



V'L;m'(^e) 



(B2) 



m 



fOO 

/o 



(B3) 



/,ph 



V'L:m'(^e) Jm'(fcsfcep) 



z/CsZ pikp 



z/Cs^: pikp 



and Jm are the Bessel functions of the first kind. Sim- 
ilarly, we follow the above procedure in cylindrical co- 
ordinates for spin waves. In that case we move to the 
reference frame used in the derivation of Eq. (Hi}, in 
which the eigenmodes are shifted relatively to the ones 
of Eq. dSD via i/';f„(K) = V'f4(K-|-kp-fc,e^). As there are 
no constraints on K, we obtain a general expression for 
the (/to) eigenmode in position space 

■Y^e^^^'^tlMP^^) (B4) 

m' 

with 

Jo 

(B5) 

and i't-^^, {K,, z) = jr^dK^ e'^^^V^f™;™' i^s, K.). 

In the case of photons, we fixed |k| of ip^^{\i) mode 
and applied the paraxial approximation in the position 
space. Hence, the z dependence of the corresponding 
cylindrically symmetric V'^mm' (z"' ^) nmst be exactly the 
one of the Laguerre-Gauss mode functions, that is, 



d_ 

dz 



pdp lg?hp,z) vf:.„.(p,^) 



/,Ph 



= 0. 



(B6) 

We verify the above claim by substituting for the pho- 
tonic modes of Eq. (|B3[) into Eq. (|B6I) with the Laguerre- 
Gaussian modes defined as 



pimip 



l\ 



w V TT y (m + l)\ 

-,2 



p 



LG^r) = e™"^ LGfV, z) 

where the reduced radial and propagation distances re- 
spectively read 



w^l + ^^,) 



LHPs), (B7) 



Ps = 



V2p 



'Vi+es 



and 



6 = 



2z 
ksW^ 



The LG modes form a complete orthonormal basis in 
the radial plane at any z, that is. 



V: 



Jdp LCrirT LG;^''(r) 



Sl.l' Sr, 



Although, in the case of atoms, neither the magnitude 
of momentum is fixed nor the paraxial approximation 
holds, it is also correct to use the Laguerre-Gaussian 
modes as the basis, but independently for each z. 

Finally, we are able to write down a general basis 
expansion of the eigenmodes of Eq. (1^^ : 



I'm' 

^fl(r) = ^cit^,,„,(z)/,f„,(r). 



?a) 
(B8b) 



with the photonic and atomic basis functions defined as 



fV^lr.\ _ Jk 



/nr) = e-=^LGr(r) /r» = 



^e*''""" LGr(rf. 



Having chosen a functional basis, we decompose the 
left-hand side of Eq. (|24p and obtain its matrix repre- 
sentation: 



dpdp'fr^^irfT -^sinc(^)Mk,K) 



I'm' 



dp dp' LGr(r)*e-*'==" I dk dK e^'^'-e'^'-' '^^^ ^'^ 



X /i(k,K)e 



iksZ 



--iku-r' 



•■"LGI7(r'), 



27rfc2 
(B9) 



fixing again the Stokes photons momentum at |k| = kg 
and substituting for the coupling function according to 
Eq. dUl), that is, 



h(k 



with h{p, z) = ^ exp 



<y^{i+ii) 



Z{z). 



The integral (jB9[) corresponds to the matrix element 
him-i'm' obtained by projecting onto the ff^ and //t^, 
basis functions. Since h{p, z) is cylindrically symmetric, 
the matrix elements are found to be diagonal in the 
mode number to, that is, 



Hm:l' 



where 



,(z;z') 



i^lm;!' m\Z] Z ) K' Oj^iYn' "/■/' (,-2^ j 



h\%{z') = {2^f pdp p'dp'hGTip.zf 



(Bll) 



X kgdkgj^{kskgp)e ' °2''(^ '-'hmikskgp') 
Jo 

xh{p',z')LGf{p',z'). (B12) 



In the last step of Eq. (|BTT|) . similarly to Eq. ((B2)) . 
we employed the paraxial approximation in the case of 
photons. As proven in Eq. (|B6|) . the result (|B12[) is 
independent of the photonic axial coordinate. Whence, 



not only /i|T^, is a function only of z', but also the z 
coordinate can be arbitrarily fixed. Taking z — z' , we 
obtain 



/^M'(^') = K / P'dp' LGrip',zrh{p',z')LG]?ip',z'), 
Jo 

(B13) 

where \s = 2'k /kg and we have used the identity 



V: rtdtJra{tp)Jra{tp')= ^^P ^'\ 

™ Jo P 



The integral in Eq. (|B13P can be computed analytically, 
as an overlap of a single Laguerre-Gaussian mode with 
a product of a Gaussian and a Laguerre-Gaussian func- 
tions [21|, that is, 



KA^)-'-fziz)i2xr^\hi^)h;' 



i+i' 



l-^L 



[1 + 2(1 + iOx + 4'] 1+"+'+'' V 1 + *e 



'-( 



X 2F1 



-I, —l';m + l;4x 



2 / 1 + *C 



1 + e? 



(B14) 



where x = (t^/w^ and 2F1 is the Gaussian hypergeomet- 
ric function. This result resembles the one calculated 
for the eigenmode problem in spontaneous parametric 
down-conversion [22]. Furthermore, the z dependence 



10 



can also be discretized in the experimentally valid re- 
gion, so that the overall matrix defined in Eq. (|B14p fully 
describes the system in a finite dimensional form. For 
example, in the cases discussed in the paper, by fixing 
the width parameter of the LG modes to w ~ \/ <Ja,zl^a 
in the ultra-cold regime (a), and w = ^jL/kg in the 
room-temperature regime (b) , we can reparametrize the 
problem and, rather than the relevant ranges \z\ <'iaa,z 
and \z\ < -J, discretize the propagation distances, ^s, in 
l^sl <6 and |Cs| < 1 respectively. 

Finally, the continuous integral eigenmode problem 
of Eq. ([5]) is reduced to performing the singular value 
decomposition (SVD) of the matrix: 

^l;l'{^) — 2^ Ckm CfcTO.;„ C^,^.;/^(z), (B15) 

k 

where the diagonality of him-i'm in "i (the cylindrical 
symmetry of the coupling function) fixes completely the 
azimuthal angle's mode number. Hence, the expressions 
for photons and spin waves eigenmodes of Eq. (jB8[) sim- 
plify to the form of Eq. (pSjl stated in the paper 

V't(r) = e-^'^^^ e^''-- e"™^ Vt (p, z). 

However, now. having computed the coefficients c^'^'.^'J^, 
we can reconstruct the effective eigenmodes as 

V-r^lp, z)-J2 <i'rn LG.'IXp, z) (B16a) 

I' 

^Up. ^) = E cr™;r™(^) LG[;Xp, zf. (B16b) 
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